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We address recent spin transport experiments in ultracold unitary Fermi gases. We provide a 
theoretical understanding for how the measured temperature dependence of the spin diffusivity at 
low T can disagree with the expected behavior of a Fermi liquid (FL) while the spin susceptiblity 
(following the experimental protocols) is consistent with a Fermi liquid picture. We show that the 
experimental protocols for extracting Xs are based on a FL presumption; relaxing this leads to 
consistency within (but not proof of) a pseudogap-based approach. Our tranport calculations yield 
insight into the measured strong suppression of the spin diffusion constant at lower T. 
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Recent measurements associated with mass transport 
[l[ and spin transport [2] in the ultracold Fermi gases and 
superfluids are of great interest principally because they 
provide detailed information about the excitation spec- 
trum, thereby strongly constraining microscopic theories. 
Equally important are their widespread implications for a 
host of different strongly correlated systems ranging from 
quark-gluon plasmas to the high Tc cuprates |3|. Com- 
mon to these materials are the short mean free paths 
which lead to "near-perfect fluidity" [ll "bad metallic- 
ity" [3| and now bad spin conductivityf?] . These latter 
are Fermi-gas-based materials with small spin diffusivi- 
ties Ds that approach the quantum limited value h/m as 
the temperature is lowered. However, the nature of the 
excitations in these Fermi gases is currently under debate 

Q and a controversy has emerged as to whether the 
associated normal state of these superfluids is a Fermi 
liquid or contains an excitation (pseudo)gap. Moreover, 
there are mixed features in recent spin transport experi- 
ments [2I] which are indicative of both pairing and Fermi 
liquid theory. 

The goal of this paper is to address these spin trans- 
port experiments by Sommer et al. '2] which measure 
the spin diffusivity Dg , and the spin conductivity as and 
thereby deduce the spin susceptibility Xs = <^s/Ds, as 
functions of temperature at unitarity. Our calculations, 
based on a pseudogap approach to BCS-BEC crossover 
theory [tI , have the aim of resolving apparent experimen- 
tal descrepancies between interpretations of Xs^ Ds and 
as- It should be stressed, however, that we focus on lower 
T (below the pairing onset T*) than in Ref. ^ where one 
can fully contrast Fermi liquid and pseudogap theories. 
Indeed, if one is to sort out whether the normal state of 
a unitary gas has a pseudogap or is a Fermi liquid, it is 



essential to address these spin (as well as mass) transport 
experiments using multiple theoretical frameworks. 

Our key observation is that one must differentiate be- 
tween the spin susceptibility as computed assuming a 
constant number of carriers in as and one which is calcu- 
lated assuming the number of carriers in as increases with 
T reflecting the presence of a pairing gap. For the latter 
case one recovers a consistent interpretation of both Xs 
and Ds , in contrast to the former Fermi liquid approach, 
which leads to this mixed interpretation. 

The physics of transport in the unitary gas is com- 
plicated because, if a pseudogap is present in the nor- 
mal state, one should accomodate both fermions and 
fermionic pairs. In recent papers we have included both 
types of quasi-particles in addressing the observed 
anomalously low shear viscosity [1] both above and below 
Tc, as well as magnetic and non-magnetic Bragg scatter- 
ing [9^. While previous studies of spin transport in the 
cold Fermi gases have employed, for example, Quantum 
Monte Carlo, Boltzmann transport and variational ap- 
proaches l3, Uj, we choose to use the Kubo formalism 
which more readily addresses conservation laws and sum 
rules. The underlying theory is a BCS-BEC crossover 
theory in a t- matrix approximation 7]. At T = 0, the 
system is in the BCS-Leggett ground state where there 
are only superconducting pairs (characterized by the gap 
Asc)- For temperatures < T < Tc, there is a comix- 
ture of superconducting pairs, noncondensed pairs (char- 
acterized by the gap Apg), and fermionic quasiparticle 
excitations. 

Essential to satisfying conservation laws and sum rules 
is the incorporation of Ward identities and collective 
mode physics. (The latter applies only to the mass trans- 
port or "electromagnetic" response.) These issues are 
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Figure 1: (a) The quasiparticle inverse lifetime tJ^ in units Ej. The values are inferred from RF experiments [T^ and are 
used as input in our microscopic calculations '8]. The temperature range is limited to O.STc < T < l.STc where the RF-derived 
lifetime data are available, (b) The spin diffusivity Da for a homogenous unitary Fermi gas in our microscopic theory. It is 
normalized by its value at Tc- (c) The spin susceptibility xf ^ computed assuming a temperature independent effective number 
of carriers and normalized by the value at Tc. At low T, Figs, (a), (b) and (c) can be directly compared to Figs. 2,3 and 4a 
respectively in Ref. 0. 



addressed in Refs. H, The presence of both con- 

densed and non- condensed pairs in transport leads to 
the usual Maki-Thompson (MT) and Aslamazov-Larkin 
(AL) diagrams, which are demonstrably consistent with 
gauge invariance. A different diagram sub-set necessar- 
ily appears in the spin response, as compared with the 
mass-transport (or electromagnetic response) ; in the for- 
mer the AL diagrams are not present, nor do collective 
mode (or phononic Goldstone bosons) enter [^, [3] ■ The 
fermionic self energy on which our work is based has be- 
come the literature standard for (T > Tc) pseudogap the- 

MM 



ones 



E(p, ia;„) 



-«7 - 



A2 



A? 



(1) 



where is the free quasiparticle dispersion and 7 is the 
fermionic inverse-lifetime associated with the conversion 



from fermions to non-condensed pairs. This general form 
for E in Eq([T|) has been applied in both experimental 17 
and theoretical [IS] radio frequency (RF) studies of the 
cold Fermi gases. In the weak dissipation limit where 
7 is small, there is little distinction between condensed 
and non-condensed pairs, while in the strong dissipation 
limit, this distinction is enforced. In this paper we find 
that these two approaches tend to converge for s-wave 
pairing and the lifetime parameters obtained from RF 
experiments 12|. For reasons of transparency, for the 



most part the equations we present are in the weak dis- 
sipation limit. 

We write the key equations first for the dc spin con- 
ductivity in the more general strong dissipation limit and 
second for the the spin-spin correlation function at gen- 
eral {cijio) (called (3gQ(q, cj)) but in the limit in which 
the non-condensed pair lifetime is relatively long: 
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where P+ = (p + ^_^^'^n + P = (p,ia;„), and 

w"*" = cj + iO^ . The quantity iuJn (j-^i) is a fermionic 
(bosonic) Matsubara frequency. The excitation energy is 

Ep = + A2 where A = ^ ^Ic + ^Ig- The super- 

cripts ± on and E'^ indicate the momentum argument 



p± ^. Here / is the Fermi function and G = G'(E) is the 
dressed Green's function. Note that while one can inter- 
pret Fsc as the usual Gor'kov Greens function reflecting 
superconducting order, there must also be a counterpart 
Fpg (discussed in detail elsewhere [3] which reflects non- 
condensed pairs). Interestingly, in the weak dissipation 
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limit the spin transport correlation functions depend only 
on the total pairing gap. When pg effects are dropped, 
these equations reduce to their usual BCS counterparts. 

Any consistent theory of the spin transport must be 
compared with the f-sum rule 



(4) 



where x//(q, = — ImQpQ/Tr and we have verified from 
Eq. ^ that this can be proved analytically. 

In the weak dissipation limit, it follows from 
Eq([2l) and Ref9 that the spin conductivity as = 



lim lim — Im^'"^ 
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Importantly, it is common to 



interpret the spin conductivity in terms of an effective 



carrier number 
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where Tg = I/7 which can be directly associated with 
p3— , the spin drag relaxation time. Here Vp = dEp/dp. 
The microscopic Kubo calculation presented here leads to 
an identification between the spin-drag lifetime, and 
the quasi-particle lifetime 7~^, since spin and "charge" 
or mass are carried by the same quasi-particles. This 
identification was also observed in previous Kubo calcu- 
lations of spin diffusion in Helium-S 19]. That the inter- 
conversion between fermions and bosons is the physical 
origin of the lifetime leads us to speculate that Ts is small- 
est for the unitary gas Q where the number of bosonic 
pairs and fermions are roughly comparable. This appears 
consistent with the findings in Ref. [3- The spin suscep- 
tibility is Xs^ — ~ ) • Alternatively, from the 
p ^ 

small uj and q hydrodynamics, it is seen that Xe = cts/ Ds- 
This latter approach is used in recent experiments 

We now turn to these experiments. In FigUfa), the 
values of inferred from RF data [12] are plotted in 
units Ep y. The lifetime increases for lower temper- 
atures since inter-conversion processes cease when non- 
condensed bosons disappear; in this way fermions become 
long lived. At the lowest T considered, this figure shows 
similar trends to Fig 2 in Ref.'^. Here one associates Tsd 
with . 

The microscopically computed spin diffusivity Ds is 
shown in Fig|IJb). Importantly we find the spin diffu- 
sivity is suppressed at low temperatures reflecting the 
suppression in the spin conductivity as as a result of the 
reduced number of carriers. As in experiment, Ds is the- 
oretically found to decrease and (at the lower T) this 
figure shows similar trends to Fig. 3 in Ref. S We es- 
timate Ds at Tc is 0.61/i/m, which when corrected by a 
factor of 5.3, estimated to account for trap effects [2|, is 
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Figure 2: (a) The effective carrier number (n/m)cfi for a ho- 
mogenous unitary Fermi gas in our microscopic theory nor- 
malized by the value at Tc- It vanishes at low temperature 
because of the pairing gap A. (a) Inset:The spin conductivity 
(Ts for a homogenous unitary Fermi gas in the microscopic the- 
ory normahzed by the value at Tc- (b) The spin susceptibihty 
Xf® for a homogenous unitary Fermi gas in our microscopic 



theory normalized by the Tc value. 



in reasonable agreement. As T — ^ wc cannot rule out 
an upturn in Dg, depending on Ts(T). 

The protocol used to deduce the spin susceptibility in 
Ref. is based on 



xf^ ^ aslDs = —Ts/Ds — = constant (6) 



n 
m 



We use the superscript FL to emphasize that this anal- 
ysis builds in a Fermi liquid interpretation by assuming 
that the carrier number is a constant in temperature and 
that no excitation gap is present. Qualitatively simi- 
lar, low T experimental behavior is shown in Fig. 4a of 
Ref. S Nevertheless, this agreement between theory and 
experiment should not be assumed to support a Fermi- 
liquid based interpretation- but rather to establish inter- 
nal consistency. [If the behavior for as is assumed to be 
Fermi liquid like, the computed Xs will also be so]. As in 
experiment, x^^ increases for low temperature and ap- 
pears incompatible with a pairing theory since the spin 
susceptibility does not vanish for low temperatures. 

This will be contrasted with the microscopic calcu- 
lation of the spin susceptibility, based on as (Eq.®), 
within a pseudogap formulation, where {n/m)cs is 
strongly temperature dependent. Plotted in Figl2ja) is 
this effective carrier number (n/m)eff as a function of 
temperature while the inset indicates the behavior for 
as, computed from Eq.®. Both as and {n/m)eS vanish 
with decreasing temperature due to the pairing gap. 

Our calculation of the spin susceptibility = cTs/Ds 
is shown in Figl2][b). This has the expected tempera- 
ture dependence associated with a pseudogap. Our key 
observation here is that one must differentiate between 
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the spin susceptibility as computed assuming a constant 
number of carriers in Us and one which is calculated as- 
suming the number of carriers in Ug increases with T 
reflecting the presence of pairing. For the latter case one 
recovers a consistent interpretation of both Xs and I?s, 
in contrast to the former Fermi liquid approach. 

We have not, in this paper proved one way or the other 
whether there is strong evidence for a pseudogap associ- 
ated with the experimental results presented in [2| ■ What 
we have established is that the indications from the spin 
susceptiblity, which are used to support Fermi liquid the- 
ory are based on Eq. © which is associated with a Fermi 
liquid like behavior. This does not establish the absence 
of a pseudogap. Other claims for this absence, based 
on thermodynamics of a balanced unitary gas [5|, can 
be countered by noting that the same thermodynamic 
power fits can be found for a non-Fermi liquid (for ex- 
ample in BCS theory below Tc). Interpretation of recent 
experiments @ which address the low T normal phase 
associated with a Fermi gas of arbitrary imbalance also 
fall in this category. The only theoretical phase diagram 
[2l| (of which we aware)- addressing where pseudogap 
and Fermi liquid phases are stable- predicts this low T 
phase should be Fermi liquid, as observed 0. 

A suggestion for establishing how to rule in or rule 
out a pseudogap in the normal phase of the unitary gas 
is discussed in Ref. |22. Essential is that one first con- 
firm the known characteristics of the superfluid phase, 
such as a suppressed spin susceptibility or entropy, and 
then establish that these features vary rather smoothly 
persisting somewhat above the transition. These obser- 
vations are the counterpart of those first applied to the 
cuprates and recent cold gas experiments [23| have sug- 
gested new techniques for measuring these potential spin 
susceptibility suppressions. 

Of great interest is the relation between spin and mass 
transport. Because spin and "charge" or mass are car- 
ried by the same quasi-particles, even in the presence of 
a pseudogap [§], spin and mass transport behave sim- 
ilarly. The analysis applied in this paper was used to 
anticipate that the anomalously low shear viscosity Q 
of the normal state should persist down to T k, 0, as 
now observed [ij. As in Ref. 24 we find the excitation 
gap is responsible for this behavior. It would similarly 
explain bad metallicity Q in the pseudogapped high Tc 
superconductors . It is striking that experiments from 
the high Tc community [isj have now tended to focus on 
the temperature dependence of the efii'ective number of 
carriers (in the presence of a pseudogap) and note that it 
will affect transport "because rieff may be changing with 
T." These commonalities highlight the importance of the 
ultracold gases as a powerful simulation tool for a wide 
class of condensed matter systems. 

In summary, in this paper we presented a theory of spin 
transport in a non Fermi liquid (FL) scenario and have 
shown that these experiments do not provide evidence 



against a pseudogap. Following the experimental pro- 
tocols for extracting the spin susceptibility we find that 
Xs appears FL like, but we emphasize this refiects a FL- 
starting point (presuming T-independent carrier number 
in as) in the protocols. As in experiment, at lower T we 
find that Ds is suppressed. Moreover Us is even more 
so, as both reflect the suppressed carrier number due to 
pseudogap effects. 
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